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ABSTRACT 



Gamma-rays from pulsars can be efficiently attenuated in 
their magnetospheres via the mechanism of single-photon 
pair production and also the exotic QED process of pho- 
ton splitting, which become prolific in fields approach- 
ing the quantum critical value of B cl - — 4.41 X 10 13 
Gauss. Recently we have published results of our mod- 
elling of strongly-magnetized 7-ray pulsars, which focused 
on the escape or attenuation of photons emitted near 
the pole at the neutron star surface in dipole fields, in 
a Schwarzschild metric. We found that pair production 
and splitting totally inhibit emission above around 10-30 
MeV in PSR1509-58, whose surface field is inferred to be 
as high as 0.7B CI . Our model pulsar spectra are consis- 
tent with the EGRET upper limits for PSR1509-58 for a 
wide range of polar cap sizes. Here we review the princi- 
pal predictions of our attenuation analysis, and identify 
how its powerful observational diagnostic capabilities re- 
late to current and future gamma-ray experiments. Diag- 
nostics include the energy of the gamma-ray turnover and 
the spectral polarization, which constrain the estimated 
polar cap size and field strength, and can determine the 
relative strength of splitting and pair creation. 
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1. INTRODUCTION 



Magnetic one-photon pair production, 7 — ► e + e~ , has 
traditionally been the only gamma-ray attenuation mech- 
anism assumed to operate in polar cap models for radio 
(e.g. Sturrock, 1971) and gamma-ray pulsars (Daugherty 
& Harding 1982, 1996; Sturner & Dermer 1994). Such 
an interaction can be prolific at pulsar field strengths, 
specifically when the photons move at a substantial an- 
gle #itB to the local magnetic field. Pair creation has a 
threshold of 2m e c 2 w 1.02 MeV for 6> kB = 90°. The 
exotic higher-order QED process of the splitting of pho- 
tons in two, 7 — ► 77 , will also operate in the high field 
regions near pulsar polar caps and until very recently, 
has not been included in polar cap model calculations. 
Magnetic photon splitting has recently become of inter- 
est in neutron star models of soft gamma repeaters (Bar- 
ing 1995, Baring and Harding 1995, Thompson & Dun- 



can 1995), mainly because of their purportedly extreme 
fields ( iS 10 14 Gauss). Splitting becomes more effective 
in competition with pair creation as a photon attenuation 
mechanism at higher field strengths (Baring 1991). 

The key property of photon splitting that renders it rele- 
vant to neutron star environs is that it has no threshold, 
and can therefore attenuate photons below the thresh- 
old for pair production, 7 — ► e + e~ . Hence the impor- 
tance of photon splitting in gamma-ray pulsar models 
clearly needs to be assessed because (i) it is an addi- 
tional attenuation process for gamma-ray photons that 
can produce cutoffs in the spectrum, and (ii) when it is 
comparable to pair production, it will diminish the pro- 
duction of secondary electrons and positrons in pair cas- 
cades while effectively softening the emission spectrum. 
Such "quenching" of pair creation can potentially pro- 
vide a pulsar "death-line" at high field strengths. About 
a dozen radio pulsars have magnetic fields, determined 
from dipole spin-down, above 10 13 Gauss. This group 
includes PSR1509-58, the gamma-ray pulsar having the 
lowest high-energy spectral turnover of ~ 1 MeV (Matz 
et al. 1994, Wilson et al. 1993, Bennett et al. 1993). Lit- 
tle attention was paid to 7 — ► 77 in pulsar contexts 
prior to the launch of the Compton Gamma-Ray Ob- 
servatory (CGRO) in 1991 because until then, the three 
known gamma-ray pulsars had estimated field strengths 
of less than a few times 10 12 Gauss. The detection of 
PSR1509-58 by the OSSE and Comptel experiments on 
CGRO provided the impetus to focus on high-field neu- 
tron star systems. 

Recent investigations of the role of 7 — ► 77 in gamma- 
ray pulsar polar cap models have been performed by 
Harding, Baring and Gonthier (1996, 1997) and Chang, 
Chen and Ho (1996). In this paper, we outline the impor- 
tance of photon splitting for gamma-ray pulsar models, 
focusing on the major features presented in Harding, Bar- 
ing and Gonthier (1997), and identifying powerful obser- 
vational diagnostic capabilities of our attenuation model 
and their relationship to current and future gamma-ray 
experiments. Principal diagnostics include the maximum 
energy of gamma-ray emission, and the spectral shape 
and the polarization below this maximum. These con- 
strain the estimated polar cap size and/or field strength 
(i.e. location of emission on or above the stellar surface), 
and can further determine the relative strength of split- 
ting and pair creation, which becomes salient for some 
subtle physics issues pertaining to splitting. As such, 
phase-resolved spectral measurements in the soft gamma- 



ray band can provide a wealth of diagnostic information, 
a goal that is readily achievable by the Integral mission. 
Future instrumentation with polarization sensitivity will 
further enhance our understanding of gamma-ray pulsars. 



2. SPECTRAL ATTENUATION IN 
GAMMA-RAY PULSARS 



Before discussing the attenuation of gamma-rays in pulsar 
magnetospheres, it is instructive to briefly review photon 
splitting. Since pair creation is widely invoked in pulsar 
models, it is appropriate here to omit a detailed discussion 
of its properties, referring the reader to Daugherty and 
Harding (1983). While 7 — > 77 is forbidden in field free 
regions, it becomes quite probable in neutron star fields, 
where B becomes a significant fraction of the quantum 
critical field B CI = m 2 c 3 /eh = 4.413 x 10" Gauss. 
Splitting is polarization-dependent in the birefringent, 
magnetized vacuum, implying that polarized photons 
emerge from an emission region. The three polarization 
modes permitted by CP-invariance in QED are _L — ►|||| , 
_L — >^_L and || — >_L|| . However, Adler (1971; see also 
Usov and Shabad 1983) devised additional restrictions, 
called polarization selection rules, by demanding absolute 
four-momentum conservation and solving the dispersion 
relations for photons in the polarized vacuum in the limit 
of weak dispersion. This limits splitting to just one polar- 
ization mode ( _L — >|| || ) below pair creation threshold and 
for B ^ B cr . Such selection rules are well-defined only 
in cases of very weak linear dispersion, and could well 
be modified in regimes of moderate or strong dispersion 
( B iJ -Bcr), or by non- linear dispersion and field non- 
uniformity effects; these modifications are not yet fully 
understood. 

The study of the physics of 7 — ► 77 has at times had 
a tumultuous history (Baring 1991, Harding, Baring and 
Gonthier 1997); the first reliable calculations of its rate 
were performed in the early 70s (e.g. Adler 1971, Papa- 
nyan and Ritus 1972), and splitting is still the subject 
of some controversy. For photon energies emc 2 with 
£ <C 1 , and fields B <C B CI , the splitting rate (e.g. 
Adler 1971) averaged over photon polarizations (e.g. Pa- 
panyan and Ritus 1972; Baring 1991) can be expressed 
as an attenuation coefficient T sp , which is the rate of 
splitting divided by c : 

T ^~ a f f 19 V mc ~5 f B sin 9 kB \ 6 m 

M£J ~ 10lT 2 J k \ B cr J ' [l) 

where ctf ~ 1/137, and #kB is the angle between the 
photon momentum and the magnetic field vectors; e is 
in units of mc 2 . Note that T sp (e) — T sp (e)R is the 
optical depth for an emission region of size R , when B 
is spatially uniform. Reducing #kB or B dramatically 
increases the photon energy required for splitting to op- 
erate in a neutron star environment, a property that also 
holds for pair creation. High field ( B <; 0.3i? cr ) correc- 
tions to the above formula (e.g. see Harding, Baring and 
Gonthier 1997, hereafter HBG97) for splitting diminish 
its dependence on B , causing the attenuation coefficient 
to saturate above B ~ 4£? cr . 

In this paper, the attenuation of photons is determined 
as in HBG97 by following photon paths in the dipole ge- 
ometry of a neutron star field, determining where they 
split or create pairs. The photons originate at different 
magnetic colatitudes 9 on the neutron star surface and 



propagate outwards, initially more-or-less parallel to the 
magnetic field. We choose the neutron star radius to be 
10 cm, and a neutron star mass of IAMq . Opting for 
points of emission on the stellar surface maximizes the 
average field strength along a photon path, thereby pro- 
ducing the highest possible optical depths for 7 — ► 77 
and 7 — > e + e~ . Emission from above the stellar surface, 
as would occur for processes such as curvature radiation 
(e.g. Daugherty and Harding 1994) or magnetic Comp- 
ton upscattering (e.g. Sturner and Dermer 1994), pushes 
the threshold for spectral opacity (i.e. the turnover) up 
in energy, since the rates for splitting and pair creation 
are increasing functions of B and energy. An important 
development in HBG97 was the inclusion of the general 
rclativistic effects of curved spacetime in a Schwarzschild 
metric, following the treatment of Gonthier and Hard- 
ing (1994). These included curved photon trajectories, 
affecting the angles photons make to the field, the grav- 
itational redshift of photon energy as a function of dis- 
tance above the neutron star surface, and an increase in 
the dipole field strength (by about a factor of 1.4 at the 
pole) above flat spacetime values. These effects all act 
to increase splitting and pair creation optical depths and 
lower the maximum energies for radiation transparency 
in the magnetosphere, typically by a factor of 2-3 be- 
low flat spacetime values. Kerr metrics were not con- 
sidered since gamma-ray pulsar periods are much longer 
than their light crossing times. 



2.1. Photon Escape Energies 

If a photon is attenuated via either absorption process 
after a distance L along its curved trajectory away 
from the pulsar surface, then we call L its attenuation 
length. Clearly from the behaviour in Eq. ([j]), and also 
for 7 — ► e + e~~ , attenuation lengths will be decreasing 
functions of the photon energy e ; in fact, at high ener- 
gies they vary as e -5 / 7 for splitting and e _1 for pair cre- 
ation (HBG97). Generally, they will be reached only after 
a photon has propagated a sufficient distance to achieve a 
significant angle to the field. This criterion is easier to ac- 
complish away from the magnetic pole due to greater field 
curvature. Hence attenuation lengths are expected to be 
decreasing functions of the magnetic colatitude 9 ; this is 
borne out in the detailed calculations of HBG97. Since 
the rates of the two attenuation processes are strongly 
increasing functions of energy, L must approach infin- 
ity at some finite energy for each process, below which 
photons are free to escape the magnetosphere. Such en- 
ergies are called the escape energies, and approximately 
delineate the energy at which spectral turnovers are an- 
ticipated; they always exist due to the r~ 3 decay of the 
dipole field. 

Fig. 1 illustrates escape energies, as determined in 
HBG97, for photons initially propagating parallel to field 
lines (#kB,o = 0°) at the stellar surface. These are 
strongly decreasing functions of B (roughly as £? -6 / 5 
for splitting and B^ 1 for pair creation when B <C B cr ) 
and 9 (~ #~ 6 / 5 for 7 — > 77 and ~ 9~ 1 for 7 — > 
e + e~ ). These dependences and their sensitivity are nat- 
urally expected to be borne out in spectral turnovers from 
a population of sources; we argue below the value of this 
strong diagnostic. HBG97 showed that if the photons are 
permitted to move at some small angle (at least around 
0.57° ) to the field initially, as might be the case for 
resonant Compton upscattering polar cap models (e.g. 
Sturner and Dermer 1994, where particles with Lorentz 



factors 7 <; 100 emit the photons), the sensitivity to 
colatitude is all but obliterated for 9^2°. Also evident 
in Fig. 1 is that for low fields (below ~ 0.3-B cr ), pair 
production escape energies are below those for splitting, 
but in high fields, splitting escape energies are lower at all 
9. Hence, one expects photon splitting is irrelevant to the 
consideration of the Crab and Vela pulsars, whose spin- 
down fields are ~4x 10 12 Gauss, but would be very im- 
portant for PSR1509-58, which has B~3x 10 13 Gauss. 
Observe that the pair production escape energies are 
bounded below by the pair threshold 2mc 2 / sin#kB, but 
photon splitting can attenuate photons well below pair 
threshold. 



Strong polarization of the continuum results from unpo- 
larized injection; clearly a polarized injection could be 
either further polarized, or depolarized by the attenua- 
tion. The spectra are determined for injection of photons 
at a small angle to the surface field. This introduces a de- 
pendence of the energy of the cutoff on the polarization 
state of the photon, which arises because of the differ- 
ent cutoff energies for photon splitting and pair creation. 
Such a property, which can be a powerful observational 
diagnostic, all but disappears when the photons start off 
parallel to the field. Note, however, that polarization 
of the continuum below the cutoff is still present when 
#kB,o = 0° . 
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Figure 1: The energy (in units of mc 2 ), below which 
photons escape the magneto sphere without splitting (solid 
curves) compared to the escape energies for one-photon 
pair production (dashed curves) as a function of magnetic 
colatitude 9 of the emission point on the neutron star sur- 
face, for different surface dipole magnetic field strengths 
(the escape energy drops as B increases). The curves, 
which are for averages over the photon polarizations, di- 
verge near 9 = because the photons are almost parallel 
to the field lines throughout their path. The dots depict 
the escape energies for the polarization mode _L— >|||| . 



2.2. Spectral Attenuation in PSR1509-58 



Here, the depiction of spectral attenuation caused by 
splitting and pair creation will be focused on the case 
of the high field gamma-ray pulsar PSR1509-58. This is 
because the "GeV" cutoffs and reprocessing in the Crab 
and Vela spectra are well-studied (e.g. Daugherty and 
Harding 1982, 1996), and as mentioned above, photon 
splitting is unimportant in these two sources. We adopt 
the high revised spin-down estimate B = 0.7B cr for 
PSR1509-58 following Usov and Melrose (1995). Photons 
are injected at the stellar surface with a canonical unpo- 
larized power-law continuum that is chosen to be con- 
sistent with the OSSE data points (e.g. see Matz et al. 
1993), and we determine emergent spectra ignoring (for 
simplicity) any photon generation by created pairs. De- 
tails of pair cascading will be the subject of future work. 
Fig. 2 shows the differential energy spectrum obtained in 
the case where both photon polarizations produce pairs, 
but only one splits ( _L— >|| || ) according to Adler's (1971) 
selection rules, so that no cascading of photons ensues. 
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Figure 2: Polarized spectra (see HBG97) for partial pho- 
ton splitting cascades, assuming unpolarized power-law 
emission (of index a = 1.6 , the OSSE best-fit value) not 
parallel to the magnetic field (9^B,a = 0-57°), at different 
magnetic colatitudes, 9, as labelled. Here only photons of 
polarization _L split, while those of either polarization 
produce pairs. The normalization of the spectrum is ar- 
bitrary. 
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Figure 3: The polarization-averaged spectrum, multiplied 
by energy squared, for different colatitudes 9 of surface 
emission, with only one mode of splitting permitted as 
in Fig. 2. The data points were obtained by the Ginga, 
Einstein and OSSE instruments, while upper bounds at 
higher energies belong to Comptel and EGRET viewings 
(see HBG97 for references). 



Fig. 3 depicts this spectrum in a "v-F v " representa- 
tion, along with observational data and upper limits. A 
wide range ( 2° Si 9 Si 25° ) of polar cap sizes match the 
observations for the chosen field strength. Restricting 
splitting to just one polarization mode ( _L — >|||| ) below 
pair creation threshold may be relaxed if modifications 
to Adler's (1971) polarization selection rules prove sig- 
nificant, as discussed above. If this arises (e.g. through 
microscopic momentum non-conservation), and permits 
all three modes of splitting to proceed, a full cascade 
ensues with several generations of splitting and polar- 
ization state switching, as described by Baring (1995). 
This pushes the photon spectrum to lower energies since 
photon splitting tends to dominate pair creation at such 
high field strengths. Large spectral peaks will result, with 
strong polarization signatures, as is illustrated in HBG97, 
and the phase space for polar cap sizes permitted by the 
observational data diminishes. The concurrent enhanced 
quenching of pair creation would then be very important 
for radio and gamma-ray pulsar models. 



3. DISCUSSION 



In the light of the attenuation results presented here, 
the diagnostic capabilites of future experiments appear 
great, and will significantly impact our understanding of 
gamma-ray pulsars. In observations of individual sources, 
if polarization measurements are not possible, it will be 
difficult to pin down the field strength and polar cap size 
separately; the energy of the cutoff clearly depends on 
both these quantities. However, if polarization detections 
yield information such as depicted in Fig. 2, where it be- 
comes possible to observationally determine the ratio of 
escape energies for splitting and pair creation, then B 
and the cap size for emission at any given height above 
the stellar surface can both be constrained at the same 
time. The spectral shape below the cutoff is a strong 
function of the number of modes of splitting that operate 
in the system. Prominent bumps should appear in source 
spectra (HBG97) unless the polarization selection rules 
that are mentioned above restrict photon splitting to just 
the _L — >|||| mode. Hence just through a determination 
of spectral shape, Integral and other gamma-ray instru- 
ments may well be able to discriminate between the op- 
eration or otherwise of the selection rules in pulsar mag- 
netospheres, thereby impacting our understanding of this 
aspect of strong field QED, a very enticing prospect. 

As a dramatic increase in the number of observed gamma- 
ray pulsars in the not too distant future is anticipated, we 
note diagnostic advances can be made with populations of 
these sources. Firstly, the polar cap scenario considered 
here makes a definitive prediction that the spectral cutoff 
energy decreases quickly with increasing spin-down field 
strength. While the Crab, Vela and PSR1509-58 provide 
a range of field strengths that suggests this behaviour, 
it is important to confirm or deny this trend with other 
sources. The recent marginal detection of PSR0656+14, 
whose spin-down field is ~ 10 13 Gauss, at around 100 
MeV (Ramanamurthy et al. 1996) conforms to this trend. 
A much larger population of detected highly-magnetized 
gamma-ray pulsars will probe this hypothesis and prob- 
ably provide the capability to discriminate between the 
polar cap and outer gap models, since the outer gap sce- 
narios (e.g. Romani and Yadigaroglu, 1995) are unlikely 
to produce this correlation with field strength. If this 
trend is confirmed, we may further have the potential 
to discern between the physical mechanisms responsible 



for the gamma-ray emission. Resonant Compton upscat- 
tering and curvature radiation produce different angular 
beaming patterns of radiation, and therefore a different 
spread of cutoff energies for given B and distributions of 
polar cap sizes. Hence distributions of key spectral pa- 
rameters can refine our modelling of these sources. Fur- 
thermore, the future full treatment of pair cascades will 
reveal the relative importance of curvature (flat) and syn- 
chrotron (steep) contributions to the gamma-ray spectra 
of strongly-magnetized pulsars. This will define a correla- 
tion between source spectral index and B , and therefore 
cutoff energy, trends that can be verified or refuted ob- 
servationally. In summary, clearly the era of the Integral 
mission will significantly advance our understanding of 
gamma-ray pulsars. 
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